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ABSTRACT 

If a weapon system is fired at a circular target and the 
impact point is distributed as a bivariate normal random 
variable, it is not particularly difficult to determine hit 
probability if the expected impact point is at target center. 
If, for some reason, the expected impact point is offset, the 
problem of determining hit probabilities becomes quite complex. 
Herein is developed a method to approximate such a hit 
probability . 
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I. 



INTRODUCTION 



The probability of hitting a specified target with a 
given weapon system is of particular interest to both mili- 
tary commanders and weapon systems analysts . Determining hit 
probability is not always a simple matter. In particular, if 
the impact points of the weapon system are distributed accord- 
ing to some probability distribution, and the expected impact 
point is not at target center, the problem of determining hit 
probabilities becomes quite complex. 

Over the past thirty years this problem has been studied 
in depth as the extensive bibliographies of F. Grubbs [Ref. 

1] , A. Eckler [Ref. 2], and W. Guenther and P. Terragno [Ref. 
3] illustrate. Several approaches to solving the problem have 
been investigated by these and other authors in an effort to 
simplify the computations required to estimate or actually 
compute the hit probabilities. 

As an example of one approach to the problem, Grubbs [Ref. 
1] uses what he calls a "single, straight-forward and rather 
simple technique" to approximate these hit probabilities in 
which he employs an approximate central chi-square distribu- 
tion with fractional degree of freedom, or a transformation to 
approximate normality. There are, of course, other methods 
available to approximate these hit probabilities, one of which 
will be discussed in the following pages. 



4 



The method which will be discussed is computationally 
easy, using simple mathematics in conjunction with a table of 
modified Bessel functions. These tables are usually avail- 
able to tlie analyst and enable him to accurately approximate 
the desired hit probabilities without the use of complex com- 
putations. It should be mentioned that tables of offset 
circle probabilities have been compiled by Rand Corporation 
[Ref. 4]. These tables give hit probabilities as a function 
of target radius and offset distance of expected impact point, 
but are not as readily available as the Bessel function tables 
used in the following procedure. It should also be noted 
that routines for Bessel functions are available for most all 
computer systems. 
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II. NATURE OF THE PROBLEM 



Consider a circular target with radius a. Place the 
origin of the coordinate system at the center of the target. 




Let be the probability density function of the 

impact point, (x,y) , of the projectile fired at the target. 
Let ( , pCj) be the mean of the impact point distribution and 
and 0*^ be the variance of the impact point distribution 

X 

in the x and y directions respectively. Assume covariance is 
zero. If the distribution of impact points follows a bivari- 
ate normal distribution, then 




I 

— e>tp 




Assume 0^ = = O" , and • Let ^ be the proba- 

bility of a hit on the target, then 
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and 



Changing to polar coordinates yields 

X ~ r cos B ^ ^ 
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and 



Jcos"ea0= 



If n is odd, then JcoS^ Pel 0 -Q 
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If n is even, then 



. Thus , 
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- — . - . Also, for n even and equal to 2m, 
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The Bessel function, denoted by Jp 



£ 

V^-i O 



VN\ (.vq+'p') ! 



, is equal to 
, where x is the argument of the 



Bessel function and p is the order of the Bessel function. 

The modified Bessel function of order zero is denoted by 

T / • \ > tlY'^ 

and IS equal to c 2^ ^ / s ^ ^ Vj.) ^ . 

2?r 

Thus . 
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so 



that 0 = 2-‘>r Ipf^) 



It was shown previously that the probability of hitting the 

target, , was equal to |d^cAv~, 

^7t<T o ^ o 
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It was also shown that 

suit, P = fr eTfe'^ itc r„ ( ) dv* 

^ ^ z-ho-"* ^ 



As a res 



or 



p = 



Q CC 

It xs only necessary to know the ratios, - and - , to 

cr (T 

determine the probability of hitting the target. 

The Rand Corporation has constructed tables giving these 



(X 



hit probabilities as a function of the two ratios, *5 and - , 

tr cr 

[Ref. 3]. 

Gilliland, [Ref. 4], shows that 



If; 

Luke, [Ref. 5], defines the function 
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Thus , 
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As a result. 






Luke, [Ref. 5], gives the expansion of in a 

series of Bessel functions . 
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These latter two equations yield the true hit probability 
for a weapon system fired at a circular target with radius a, 
and normally distributed impact points with expected impact 
point at and . 

A good approximation to the true hit probability can be 
had by considering only the first few terms of the infinite 
series in the appropriate equation. That is 
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R >(^ 



where rv is, hopefully, relatively small. If KN is less 
than five, then the approximation v/ould be beneficial to 
analysts and would be quite simple to use. If n is 
greater than five, the calculations required to get a reason- 
ably accurate approximation to hit probability become 
prohibitive . 
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III. CONCLUSIONS 



The accuracy of the above approximations was checked on 
an IBM 360 computer. The program used the FORTRAN language 
and two IBM subroutines, 10 and INUE, which computed the 
modified Bessel functions. This program is available through 
the Operations Analysis Department of the Naval Postgraduate 
School in Monterey, California. 

The computer program was basically an iterative process 

in which different values of the ratios ® and ^ were set 

0" CT 

and the hit probabilities computed while the number of terms 
from the appropriate infinite series varied from one to ten. 
The approximate hit probabilities thus obtained were then 
compared to the true hit probabilities from the Rand tables 
[Ref. 4]. 

The tv/o aspects of primary interest were (1) how accurate 

are the approximating equations in the various ranges of ^ 

and 2* , and (2) how many terms of the infinite series in the 

or 

approximating equations must be considered before a reasonable 
degree of accuracy is obtained. The specific values of the 

ratios 3 and ^ ranged from .5 to 5.0 in the computer program. 

«T <T 

For the equation for R<a, the ratio B was initialized at the 

0 " 

value .5 and the ratio ^ was allowed to vary between .5 and 

5.0. The value for B was then incremented to 1.0 and the 

<T 

process repeated with values of ^ greater than or equal to 

1.0. This procedure was continued until both B and were 

O" 



12 



equal to 5.0. For each value of B- and ^ the approximate 

CT O* 

hit probability was computed using only the first term of the 
infinite series initially / and then again for each additional 
term used, up to ten terms. In this manner it was possible 
to see how rapidly the approximating equations converged to 
the true hit probability. 

For the equation for R>a, the same type of procedure was 

used, holding the ratio S' constant while the ratio S varied 

CT <T 

through the appropriate values. As before, the ratio S- was 
initialized at .5 and then incremented by .5 each time the 

ratio ? had run its range of values and the hit probabilities 

O" 

computed. 

Representative values of hit probabilities obtained 
through the approximating equations are contained in Appendix 
A. Figure 1 compares approximate hit probabilities with ac- 
tual hit probabilities for R<a, and Fig. 2 does the same for 
R>a. In both cases the values shown were computed using the 
first five terms of the infinite series from the appropriate 
approximating equation. 

The convergence of the approximating equation values to 
the true hit probabilities is shown in the two examples be- 
low. Example one shows the convergence of the approximating 
equation for R<a as the number of terras from the infinite 

series increases . Values of B and ^ were selected to show 

O- O' 

both rapid convergence and relatively slow convergence. Ex- 
ample two shows convergence of the approximating equation for 
R>a. These examples along with Appendix A point out that the 
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approximation of hit probability for R<a is optimistic while 
the approximation of hit probability for R>a is conservative. 



Example 1 (R<a) 




B = 
0- 


.5, ^ = 1.0 

0" 



Approximation using; 

one term .2089828 



two terms 


.1108707 


three terms 


.1047551 


four terms 


.1045011 


five terms 


.1044937 True probability; .104491 


six terms 


.1044937 


seven terms 


.1044937 


eight terms 


.1044937 


nine terms 


.1044937 


ten terms 


.1044937 


^ _ 
cr 


4.5, ^ = 5.0 

O' 



Approximation using; 

one term .9253550 



two terms 


.8596848 


three terms 


.8044760 


four terms 


.7601165 


five terms 


.7260436 True probability; .653154 


six terms 


.7010158 



seven terms .6834279 

eight terms .6715976 



nine terms 


.6639764 


ten terms 


.6592714 
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1.0 



Example 2 (R>a) 

S- = 5 ^ = 

O' 

Approximation using: 



one term 


.0690204 


two terms 


.0732899 


three terms 


.0734669 


four terms 


.0734724 


five terms 


.0734725 


six terms 


.0734725 


seven terms 


.0734725 


eight terms 


.0734725 


nine terms 


.0734725 


ten terms 


.0734725 


~ = 4 

<r 


5 = 

' 0“ 


roximation using: 


one term 


.0656702 


two terms 


.1208790 


three terms 


.1652384 


four terms 


.1993113 


five terms 


.2243391 


six terms 


.2419270 


seven terms 


.2537574 


eight terms 


.2613786 


nine terms 


.2660835 


ten terms 


.2688691 



5.0 



True probability: .073473 



True probability; 



,272201 



Accuracy of the approximating equations can be improved 
in all cases by considering more terms of the appropriate in- 
finite series. If it is necessary to use more than five 
terms of the series, however, the mathematic manipulations 
involved become impractical to do by hand. Accuracy can also 
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be improved when - is equal to - . From Appendix A it is 

CT , CT 

apparent that when ^ is equal to ^ and both are greater 
than 3.0, accuracy is the poorest. To more accurately ap- 
proximate the true hit. probability in this case it is suffi- 
cient to compute the arithmetic mean of the values obtained 
from both equations (R<a and R>a) . For example, for 

^ = ^ = 5.0, the true probability is .460. The equation 

cr <j 

for R<a yields an approximation of .643, while the equation 
for R>a yields an approximation of .325. The arithmetic mean 
of these two values is .484, which is more accurate than 
either of the two approximations. 
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IV. SUMMARY 



If the impact point of a weapon system is distributed as 
a bivariate normal random variable, the task of computing the 
hit probability for a circular target when the expected 
impact point is somewhere other than at target center is for- 
midable . The two approximating equations that have been 
developed 



are relatively accurate approximations to true hit proba- 
bilities using only the first five terms of the series. 
Greater accuracy can be achieved by considering more terms. 

To have an accurate approximation to such a hit probability 
is beneficial to both weapon systems analysts and military 
commanders, because of the difficulties involved in comput- 
ing actual hit probabilities. The above equations yield a 
good approximation to hit probabilities using simple mathe- 
matics and a table of modified Bessel functions, which is not 
difficult to obtain. 




for R, 




for 
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APPENDIX A 



COMPARISON OF HIT PROBABILITIES 



or 


R 1 . 

app 


0 

act 


2 . 

app 


0 

act 


3 . 

app 


0 

act 


4 . 

app 


0 

act 


5 . 

app 


0 

act 


1.0 


.267 


.267 


* 


* 


* 


* 


* 


* 


* 




O 

• 

CM 


.731 


.731 


.410 


. 39 7 


* 




* 


* 


* 


* 


3.0 


.956 


.956 


.788 


. 786 . 


498 


.433 


* 


* 




* 


4.0 


.997 


.997 


.966 


.966 . 


814 


.803 . 


578 


• 450 


* 




5.0 


.999 


.999 


,998 


.998 . 


970 


.969 . 


8 35 


.813 


.643 


.460 






Figure 1 . 


Tab le 


for 


R Less 


Than a 








a. 

^ app 


0 

act 


2 . 

app 


0 

act 


3 . 

app 


0 

act 


4 . 

app 


0 

act 


5 . 

app 


0 

act 


(T 

1.0 


.267 


.267 


* 


«■ 


* 




* 






* 


2.0 


.082 


.082 


.39 3 


.39 7 


* 


«■ 


* 


«■ 


■jt 




3.0 


.011 


.011 


.113 


.113 . 


399 


.433 


* 


«• 


* 


* 


4.0 


.001 


.001 


.015 


.015 . 


121 


.126 . 


366 


.450 


* 


* 


5.0 


.000 


.000 


.001 


.001 . 


016 


.017 . 


12 0 


.133 


.325 


.460 






Figure 2. 


Table 


for 


R Greater 


Than 


a 





These tables compare the approximation to the actual hit 
probabilities (app) and the actual hit probabilities (act) 
themselves. Only the first five terms of the infinite series 
described on page n were used for the approximations. Accu- 
racy can be improved in all cases by considering more terms 
of the series. 
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